Mon. Not. R. Astron. Soc. 000,11]-?? (2004) Printed 2 February 2008 (MN WT^ style file v2.2) 



An XMM-Newton observation of the nova-like variable UX 
UMa: spatially and spectrally resolved two-component 
X-ray emission 

G.W. K. Mukai^f, B.J.M. Hassall^ T. Naylor^ and J.H. Wood^ 

^ CEA/Saclay, Service d'Astrophysique, L'Orme des Merisiers, Bat. 709, 91191 Gif-Sur-Yvette Cedex, France 

^ MPE Garching, Giessenbachfie, 85478 Garching, Germany 

^Code 662, NASA/Goddard Space Flight Center, Greenbelt, MD 20771, USA 
, ^Centre for Astrophysics, University of Central Lancashire, Preston PRl 2HE 
' ^School of Physics, University of Exeter, Stocker Road, Exeter EX4 4QL 
' ^ Keele University, School of Chemistry and Physics, Keele, Staffordshire ST5 5BG 



Submitted 2003 October 22 



ABSTRACT 

In the optical and ultraviolet regions of the electromagnetic spectrum, UX Ur- 
sae Majoris is a deeply eclipsing cataclysmic variable. However, no soft X-ray eclipse 
was detected in ROSAT observations. We have obtained a 38 ksec XMM-Newton ob- 
servation to further constrain the origin of the X-rays. The combination of spectral 
and timing information allows us to identify two components in the X-ray emission 
of the system. The soft component, dominant below photon energies of 2 keV, can be 
fitted with a multi-temperature plasma model and is uneclipsed. The hard component, 
dominant above 3 keV, can be fitted with a kT ^ 5 keV plasma model and appears to 
be deeply eclipsed. We suggest that the most likely source of the hard X-ray emission 
in UX UMa, and other systems in high mass transfer states, is the boundary layer. 

Key words: Stars: binaries: eclipsing — stars: novae, cataclysmic variables — stars: 
individual (UX UMa) — X-rays: stars. 



1 INTRODUCTION 

In non-magnetic cataclysmic variables (CVs), a white dwarf 
primary star is accreting material, via an accretion disc, 
from a low mass quasi-main sequence secondary star. The 
disc instability theory predicts a limit cycle behavior when 
the mass input rate from the secondary is low: mass gradu- 
ally accumulates in the disc at low optical brightness (qui- 
escence), then rapidly accretes onto the white dwarf during 
outburst when the disc is hot and br ight in the op tical and 
in the UV (see a recent review bv lLasotal boOll) . This is 
widely accepted as the underlying cause of the dwarf nova 
phenomenon that is seen in the majority of CVs. In contrast, 
the mass input rate may be high enough in some CVs for 
the accretion disc to achieve a steady state. This is thought 
to be the explanation for nova-like systems, which resemble 
dwarf novae in outburst at all times. 
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Simple models suggest that the boundary layer between 
the primary and the accretion disc should be extremely 
luminous, emitting up to half the accretio n energy. Sinc e 
most of this is expected to be in X-rays JPringldll977l) . 
X-ray observations have become an essential tool for the 
study of accretion discs in CVs. It was discovered early 
on that dwarf novae in quiescence were moderately bright 
hard (2-10 keV) X-ray sources. In outburst, however, the 
hard component generally became weaker, and sometimes 
bright soft (<1 keV) emission was observed. Combining ex- 
tensive Ein stein X-ray data on CVs with theoretical con- 
siderations, iPatterson fc Ravmondl (Il985al lbl) linked this to 
the optical depth in the boundary layer. For low accretion 
rates (appropriate for dwarf novae in quiescence), the bound- 
ary layer is hot (kT~10 keV) and optically thin, cooling 
via hard X-ray bremsstrahlung emission, while for high ac- 
cretion rates (dwarf novae in outburst and nova-like sys- 
tems), it becomes optical thick, emitting a blackbody-like 
spectrum with kT~10-50 eV. Optically thin, outer regions 
of the boundary layer can explain the residual hard X-ray 
emission. The recent X-ray, EUV . and optical observat ions 
of SS Cygni through an outburst iWheatlev et alJl2003l) are 
a recent confirmation of this basic picture. 
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Table 1. Log of XMM-Newton observation 0084190201 of 
UX UMa. 



Instrument 


Start 


End 


Data Mode 


MOS 1&2 


01:07 


14:33 


Thin Filter, Full Frame 


pn 


01:30 


14:23 


Thin Filter, Full Frame 


RGS 1&2 


01:06 


14:40 


Spectro+Q 


OM 


01:37 


14:10 


UV Grism 1 



All times in TT on 2002 June 12. 

Timing observations of eclipsing systems are invaluable 
as a method for determining the source of the X-ray emis- 
sion. Observations of high inclination dwarf novae in qui- 
escence have shown evidence for a co mpact X-ray source 



that is eclipsed with the white dwarf (IWood^^Ll 
Teeselin3 1 1997L iMukai et al]ll997L I Pratt et al.l 



n 



1995a 



Ramsav et alj200iri . although the X-ray emission in OY Car 
may be from a polar region rather than from an equa- 
torial boundary layer ( Whoatley & West 2003). However, 
the soft X-ray emission from the dwarf nov a OY Car in 
super out burst ^ was found to be uneclipsed iNavlor et alJ 
|l983)- Based on a comparison with lower inclination (non- 
eclipsing) systems, they advocated a two-component ori- 
gin for the soft X-ray emission from dwarf novae in out- 
bursts: a compact region, presumably the boundary layer, 
and an extended region, possibly due to the scattering of 
emiss ion from the former in an accretion disc corona or a 
wind. lNavlor et al] (^8^ argued that the white dwarf pho- 
tosphere and the boundary layer are hidden from our view 
at all orbital phases in OY Car in superout burst, due to the 
thickened rim of the accretion disc. Later observations of 
OY Car in superoutbursts confirm the lack of eclipse and 
strengthen the scattering i nterpretation jPratt et alJll999bl : 
iMauche fc RavmondllioOll) . 

XMM-Newton provides several key advantages in the 
study of eclipsing CVs, such as the high effective area of 
the EPIC instruments both at low and high energies. We 
have chosen the nova-like system, UX UMa, as our target 
for this study, which avoids the logistical complication of 
scheduling an observation at short notice during an out- 
burst of a dwarf nova. UX UMa is the brightest eclipsing 
nova-like system in the optical and in the X-rays, and has 
been well studied in ultra violet in recent years. From UV 
data. lBaptista et alj lll995l) claimed the detection of a white 
dwarf eclipse, which howeve r has been brought into ques- 
tion by'Fronina- et al.' (200^. A ROSAT PSPC observation 
bv lWood et al.. (.1995b.) showed no soft X-ray eclipse. 



2 OBSERVATION 

UX UMa was observed with XMM-Newton on 2002 June 
12 for approximately 45 ks (see Table0 for details). Peri- 
ods of high background in the EPIC data due to soft pro- 
to n solar flares wer e remo ved using the method described 
in iPratt fc ArnaudI i2003ri . There is one large flare in the 
middle of the observation, otherwise the data are relatively 



clean^. Useful exposure times after cleaning were 37895s, 
38412s and 28736s for MOSl, M0S2 and pn, respectively. 

For the MOS cameras, PATTERNS 0-12 were selected; 
we use only the well-calibrated single events (PATTERN 0) 
in the pn analysis. In addition, we select FLAG=0 to exclude 
events at CCD edges and around bad pixels. 

EPIC source light curves were extracted in Is bins from 
a circular region of radius 30". The background light curves 
were accumulated in Is bins from an annular region centred 
on UX UMa, with inner and outer radii of l'.45 and 5'. 3, 
respectively. Other sources were identified by eye and re- 
moved before the background products were accumulated. 
The EPIC spectra were extracted in the same regions as 
described above. Response and effective area files were pro- 
duced using the SAS tasks rmfgen and arfgen. 

For the RGS data, the standard pipeline-produced 
spectra and response matrix files were used. The 'fluxed' 
RGS1-I-RGS2 spectrum shown in Figure|5| was produced 
from these files using the SAS task rgsf luxer. 



3 RESULTS 

3.1 EPIC spectral characteristics 

The EPIC spectra of UX UMa are shown in FigureEl It is 
clear that the underlying emission is line-rich: the Fe Ka line 
is clearly visible, so is the "bump" around 1 keV indicative of 
the unresolved Fe L complex and other lines. Since we wish 
to optimise the choice of spectral bands for the extraction 
of the light curves (Sect. l^?^ . we will characterise the orbit- 
averaged spectrum of the system with the simplest plausible 
model and use this as a constraint for the choice of light 
curve bands. 

In view of the line-rich nature of the spectrum, we start 
with a single-temperature mekal model jMewe et alJll98l 
119861: iLiedahl et al.lll995l) for thermal emission from coUi- 
sionally excited plasma. This model fails to fit the data 
(x?=3.6), and predicts more fiux than observed around 2 
keV and less flux than observed in the 1 keV bump and at 
high energies. Even a multi-temper ature plasma does not 
fit well: with a cemekl (Sinaii ct alj |l996h model, in which 
the difi'erential emission measure is a power-law of plasma 
temperature, we obtain xS=l-6 for 563 degrees of freedom 
(i.e., formally unacceptable). Moreover, this fit is achieved 
with best-fit abundances of 3.9 times Solar. While there is 
no compelling reason to believe UX UMa has exactly Solar 
abundances, a factor of 3.9 seems too extreme to be believ- 
able. Additionally, systematic residuals are evident (see mid- 
dle panel of Figure^ even with the overabundant cemekl 
model. 

We have therefore investigated the soft and hard parts 
of the spectrum separately. The strength of the Fe Ka line at 
6.7 keV suggests a strong contribution from thermal plasma 
with kT in the 1-10 keV range (where Fe is mostly He-like). 
We thus experimented with the low energy cutoff when fit- 
ting a single temperature mekal model to the hard part of 
the spectrum. When the low energy cutoff Ec > 2 keV (Ta- 
blel^J, the fit becomes acceptable: including energies below 



^ For the purpose of this paper, normal outbursts and superout- 
bursts can be considered together. 



^ Unfortunately, this flare causes us to lose one of the 3 eclipses 
covered by these data. 
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Figure 1. The XMM-Newton EPIC spectrum of UX UMa, as determined from simultaneous fit to PN, MOSl, 
and MOS2 data: (Left) The observed PN spectrum plotted with the best-fit model (see text) convolved with 
instrument response. (Middle) The residuals, in the form of 5x for single temperature mekal (top), cemekl 
(middle), and our best-fit model (bottom). (Right) The two components of our best-fit model are shown 
overplotted on deconvolved PN data. Units are fiux units. 



Table 2. Results of single temperature mekal model fits to the 
hard component in the EPIC spectra. 



Table 4. Lines detected in the RGS spectra. 



Energy Range 


X^/D.O.F. 


Nh 


kT 


Norm. 




(keV) 




(cm^^) 


(kcV) 






5-10 


1.25/ 54 


1.0x10^2 


5.6 


4.2x10" 


-4 


4-10 


1.04/ 90 


6.9x10^2 


5.5 


4.8x10" 


-4 


3-10 


1.06/125 


5.6x10^2 


5.5 


4.4x10" 


-4 


2-10 


1.12/171 


3.8x10^2 


6.0 


3.8x10" 


-4 


1.5-10 


1.48/216 


1.7x10^2 


7.7 


2.8x10" 


-4 


1.0-10 


2.51/318 


0.0xl022 


9.8 


1.8x10" 


-4 



a normalization of 4.4lQ g x 10 ^ 
3.2x10"^^ ergss"^cm 



this threshold leads to a rapid deterioration in the goodness 
of the fit, and a drastic change in the parameter values. With 
Ec=3 keV, the best fit values (90% confidence level errors) 
are Nh=(5.6 ± 1.7) x 10^^ cm-^ kT=(5.5 ± 1.1) keV, with 

and a 3-10 keV flux of 
^ (Solar abundances assumed). Note 
the high value of Nh, which would indicate little contribu- 
tion of the hard component below ~2 keV. The temperature 
is constrained by the high energy cut-off and by the domi- 
nance of the 6.7 keV line. 

A single-temperature, solar-abundance mekal model 
fails to fit the spectrum of UX UMa in the soft (0.2-2 keV) 
range ixl=3-5 for 394 de grees of freedom) . Although it does 
much better (Xi5=l-5) when abundances are allowed to vary, 
the required abundances are 0.08 Solar, a value that would 
make fitting of the hard range extremely problematic. A bet- 
ter solution might be to use multiple-temperature plasma 
models. Here we have adopted cemekl, resulting in xS = l-2 



(393 degrees of freedom) with Nh=7.8I^:9 x 10^^ cm"^, 
power law index (for differential emission measure distri- 
bution) a = 0.44to :[g, kTniax=5.lti;2 keV, a normaliza- 
tion of l.Stol X 10"*, with a 0.2-2.0 keV flux of 2.3x10"" 
ergss^^cm^^ (Solar abundances assumed). 

We show the result of the combined flt of these two 
models over the 0.2-10 keV range in Figure^and summa- 
rize the parameters in Tabled We obtain x^=l-14 for 568 
degrees of freedom. The best-flt model to the soft part of 



Species 


Lab. A 

(A) 


Obs. A 

(A) 


Flux 

(10~^ photons s~-'^cm~2) 


Fe XVII 
Fe XVII 
O VIII 
O VII 


15.01 
17.10 
18.97 
21.60 


15.01±0.03 
17.09±0.04 
19.00±0.02 
21.60±0.04 


q Q+3.6 
11.2+11 



the spectrum makes a small contribution to the 3-10 keV 
band and hence the parameters of the model fit to the hard 
part of the spectrum are somewhat modified from the results 
shown in Table|5| 



3.2 RGS spectral characteristics 

We present the combined RGS spectrum of UX UMa in 
Figure|5| Continuum is weakly detected at best, while we 
detect 4 lines with some confidence. We have fitted the first- 
order spectra from both RGS 1 & 2 simultaneously, after 
rebinning the raw channels into groups of 4. We have used 
a power-law continuum model with Gaussians to obtain a 
rough characterisation of the lines (Table^J- Due to concern 
over the low signal-to-noise ratio in the continuum, and the 
large contribution of the background to the total counts, we 
have experimented with several different rebinning schemes. 
We can confirm that the best flt values remain within the 
formal error ranges quoted in Table^] These 4 lines can be 
identifled with Fe XVII 15.01 and 17.10 A, (although the Fe 
XVII line at 17.05 A cannot be excluded as the identification 
of the latter), O VIII Lya (18.97 A), and the resonance 
line of the O VII triplet (21.60 A). The intercombination 
(21.80 A) and forbidden (22.10 A) lines of O VII appear 
to be significantly weaker. Line widths are not significantly 
detected, but we can establish weak upper limits (full-width 
at half-maximum, FMHM, of 2000 - 3000 kms"^). 
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Table 3. Results of EPIC spectral fits in the 0.2-10 keV band. 



Comment Nh a kT/kT^ax Flux (10"^^ ergss-^cm'^) 

{cm-2) (keV) 0.2-2 keV 3-10 keV 0.2-10 keV 

Soft 7.8^4;^ X 10^8 0.44to ^6 5.115:2 2.28 0.37 2.87 

Hard S.lt^'j x 10^2 5.71^-2 0.0 2.85 3.05 

(Unabsorbed) ' ' 3.36 3.76 8.23 



■1.1 ,1 I 



3.0 10.0 keV 



o,o!|lfl 



Wavelength [AJ 

Figure 2. The XMM-Newton RGS spectrum of UX UMa, plot- 
ted against wavelength, in the 6-30A range. The four highest 
peaks are found at wavelengths of two strong lines of Fe XVII, 
O VIII Lyo, and the resonant line of O VII (wavelengths of the 
intercombination and forbidden lines are also indicated). 



3.3 Folded Light Curves 

We have seen in Sect.OUlhow the entire 0.2-10. keV EPIC 
spectrum can be adequately fitted with a two component 
model. Considerable uncertainties remain in the model pa- 
rameters, and the low energy part of the spectrum can also 
be described by more complex models (e.g., partial covering 
or partial ionisation) but this simple and plausible spectral 
fitting serves to motivate the choice of bandpasses for the 
light curve extraction. 

We have thus extracted the light curves of UX UMa 
from the three EPIC cameras in the energy bands 0.3- 
2.0 keV and 3.0-10 keV. We then folded each of the 6 
light curves on the o rbital period, using the ephemeris of 
iBaptista etlj] ^9^^. The folded light curves from the 
three cameras were then averaged to produce a single folded 
light curve per energy range. These are shown in Figure|21 
The light curves, particularly in the hard band, show consid- 
erable variability, often in the form of fiare-like events lasting 
several hundered seconds. These are not averaged out in the 
folded light curve because the observation covered only a 
few cycles. 

We confirm the ROSAT resuh JWood et alll995bft that 
no eclipse is seen in UX UMa in the soft X-rays (<2.0 
keV). However, there does appear to be an eclipse in the 
hard (>3 keV) band, at phase 0.0 as defin ed by optical/UV 
eclipse. According to lBaptista et all (^9^, the white dwarf 




3 iFroning et"ai] ^2003^ found it necessary to shift their 2001 FUV 
light curves earlier in phase by 0.006 cycles. We have not applied 
such a shift. 



1.0 

Orbital Phase 



Figure 3. The XMM-Newton EPIC light curves of UX UMa, 
folded on the orbital ephemeris of Baptista et al. (1995) in 50 
bins per cycle, in two energy bands. The 0.3—3 keV light, shown 
in the lower panel, is consistent with no eclipse. The 3-10 keV 
light curve (upper panel), however, suggests the existence of a 
deep eclipse. 



eclipse width (mid-ingress to mid-egress) is 0.053±0.001 cy- 
cles, while totality (2nd to 3rd contact) is a little over 0.04 
cycles. To investigate the significance of a hard X-ray eclipse 
with an assumed duration similar to that of the UV eclipse, 
we have re-folded the data in 25 bins per cycle, with one 
of the bins centred on the nomi nal phase 0.0 using the 
ephemeris of iBaotista et al] 1^9^ . The "eclipse" bin, (cor- 
responding to the two lowest points in Figurel^J has a count 
rate of 2.5±1.7 xlO~^ ctss~^ (the error from propagation 
of counting errors), while the other 24 bins have 34.4±8.5 
xlO"'' ctss"'^ (where the error is derived from the standard 
deviation of t hese bins, to take in to account source variabil- 
ity - see e.g., IWood et alJll995bft . We thus appear to have 
detected an eclipse at the 3.8 cr level. 

The hard X-ray eclipse thus appears real, but we are 
unable to place tight constraints on the parameters. We can 
reproduce the eclipse shape with a wide range of parame- 
ters (e.g., a mid-eclipse phase of 0.0-0.01, a mid-eclipse fiux 
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between 0-20% of out-of-eclipse flux, an ecUpse width of 
0.04-0.06 cycles, and eclipse transitions of less than ~0.01 
cycles) . 



4 DISCUSSION 

The combination of the spectral and timing results strongly 
suggests that there are (at least) two components in the 
X-ray emission from UX UMa. The soft component is unab- 
sorbed and uneclipsed and therefore must originate from an 
extended region. The hard component, which we have newly 
discovered, is heavily absorbed and eclipsed. 

The soft component in UX UMa probably has the same 
origin as that of OY Car in superoutburst llNavlyr et all 
Il988l: IPratt et alJ Il999 b': 'Mauchc & Ravmond"200lD, and 
perhaps that of DQ Her iMukai c t al. 2003). That is, the 
likely origin is scattering of centrally generated X-rays in an 
extended region. The prominence of the resonance compo- 
nent of the O VII triplet in t he RGS spectrum (F i gurei3 l 
supports this interpretation. iMaiiphe fc R.avmondl (120011 ) 
make a strong case that the scattering medium is the accre- 
tion disc wind in the case of OY Car in superoutburst. This 
may well be the case in UX U Ma, since the system is known 
to ha ve a strong wind (^Baptista et aljll99 5: Froning ct alJ 
l2003l) . Although we have not been able to confirm this by 
resolving the X-ray lines, our upper limits for their widths 
(FWHM of 2000 - 3000 kms"^) are consistent with the FUV 
line widths (FWHM ^1800 kms^^ Froning et al. 2003). 

The inferred luminosity of the hard component is 1.26 x 
10'^^ ergss~^ in the 0. 2-10 keV band for a distance of 345 pc 
iSaptista et al.lfl995fl . This value is similar to those found 
for the residual hard X-ray emi ssion component in non- 
eclipsing dwarf novae in outburst. iBaskill et al.l ll200lll infer 
a bolometric luminosity of ~ 6 x 10^° ergss"^ for Z Cam 
in outburst from ASCA data (where there is an additional 
softer component, which we would not see in UX UMa due 
to the absorber), while the hard X-ray component of SS Cyg 
in outburst reaches a minimum at about ^ 3 x 10'^^ ergss"^ 
in 3-20 keV in the RXTE data iWheatlev et al . 2003'). In 
all three systems, the existing spectra can be fitted with a 
single temperature thermal model with kT in the 5-10 keV 
range. These similarities both in luminosity and in temper- 
ature suggest that the hard X-ray component that we have 
discovered in UX UMa is the same as the hard X-ray emis- 
sion seen in dwarf novae in outburst. 

The detection of a hard X-ray eclipse allows us to rule 
out an extended origin for this em ission, such as was sug - 
gested (for SS Cyg in outburst) bv IWheatlev et al.l J2003f) . 
If the emission is coming from a compact region close to the 
white dwarf, the most natural source at this position would 
be the boundary layer. At first sight, this may seem an odd 
suggestion, since the boundary layer is normally thought of 
as a source of soft X-rays, and in an eclipsing system, is 
hidden behind the disc rim. However, UV observations sug- 
gest that a disc rim can supp l y a co lumn density of ~ 10^^ 
cm~^ in OY Car lHorne et alj 1^9^, and that the absorber 
in UX UMa appears to result i n deeper absorption in the 
UV than that seen in OY Car jBaotista et al.lll998l) . Fur- 
thermore, our spectral fits suggest a column of this order 
of magnitude. Pending more detailed analysis, we suggest 
that the same material can be responsible for both X-ray 



and UV absorptions. Crucially, such a column is sufficient 
to extinguish a soft X-ray source, but will not significantly 
affect X-rays above 5 keV. Thus these X-rays may well be 
coming from the more tenuou s parts of the boundary layer 
jPatterson fc Ravmondlll985afl . 

We require more sensitive hard X-ray observations of 
UX UMa (such as a longer observation with XMM-JVewton) 
to use the eclipse light curves to constrain the nature of the 
hard X-ray emission in these systems. At the same time, 
better quality spectra of the soft component are required to 
study the kinematics of the soft X-ray scattering region in 
UX UMa. 
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